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ABSTRACT 

We present the first results from our survey of stellar substructure in the outskirts of M31 using the 
Advanced Camera for Surveys (ACS) on board the Hubble Space Telescope. We discuss the stellar 
populations associated with five prominent stellar overdensities discovered during the course of our 
panoramic ground-based imaging survey with the Isaac Newton Telescope Wide-Field Camera (INT 
WFC); a sixth pointing targets a region of 'clean' halo. The colour-magnitude diagrams (CMDs), 
which contain between ^^10, 000-90, 000 stars and reach four magnitudes below the horizontal branch, 
reveal clear variations in morphology between most fields, indicating that the age and/or metallicity 
mix of stars is not constant at large radius. This directly confirms the existence of large-scale pop- 
ulation inhomogeneities within the halo of M31 and lends further support to the notion that M31 
has formed, at least in part, through satellite accretions. We find a striking similarity between the 
populations of the giant stellar stream and those of another overdensity, the NE shelf, which lies 
north-east of the galaxy center. If these overdensities are associated with the same population, then 
the difference in their red clump magnitudes implies the NE shelf lies in front of the stream by several 
tens of kpc, in good agreement with recent orbit calculations for the stream progenitor. 
Subject headings: galaxies: formation - galaxies: evolution - galaxies: structure - galaxies: halos - 
galaxies: individual (M31) - galaxies: stellar content 



1. INTRODUCTION 

The INT WFC survey of M31 has led to the discovery 
of significant stellar substructure in t he halo and outer 
disk of our nearest g i ant neighbour (llbata et ahl [22^ 
Ferguson et al.' * 200a iMcConn achie et all I2003L 12001 
Irwin et a l. 2005). At the present time, the survey con- 
sists of more than 40 square degrees of V- and i-band 
imaging to a depth sufficient to resolve stars to 3 mag- 
nitudes below the tip of the red giant branch (RGB). 

Few constraints currently exist on the nature and ori- 
gin of the stellar substructure in M31. Within the pop- 
ular ACDM model for galaxy formation, substructure is 
expected in the outer regions of galaxies as they continue 
to grow from the accretion and tidal disruption of satel- 
lite companions. Recent calculations suggest that mas- 
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sive disk galaxies, like the Milky Way and M31, could 
have cannibalised ~ 100 dwarf galaxy-like systems over 
their lifetimes (Bullock & Johnston 2004). Since these 
accreted satellites are likely to have experienced a range 
of star formation and chemical evolution histories, one 
expects their tidal debris to be distinct in terms of stel- 
lar populations. If this picture is correct, then the spatial 
substructure in M31 should be accompanied by stellar 
population inhomogeneities. 

In order to gain insight into the origin of the M31 
halo substructure, we are conducting a detailed study of 
the associated stellar populations using the HST/ACS. 
Our INT WFC survey indicated intriguing variations in 
the mean RGB colour across the halo, which could re- 
flect changes in metallicity and/or age. HST is required 
to confirm these variations and, via detection of fainter 
and more evolved populations such as the red clump 
and horizontal branch, provide constraints on their na- 
ture. In this Letter^ we present the first results from 
our HST/ACS survey which establish clear differences, 
as well as some similarities, between the populations in 
various substructures. 

2. OBSERVATIONS 

Observations at six locations in the outskirts of M31 
were obtained during Cycle 11 with the Wide Field Cam- 
era of the ACS as part of GO#9458. Our fields sample 
the five most pr ominent stellar overdensities in the maps 
of iFerguson et al.l ( 2002) as well as a region of appar- 
ently 'clean' halo (Minor Axis) situated 20 kpc along 
the southern minor axis (see Figure 1). The main char- 
acteristics of the stellar overdensities can be summarized 
as follows: 

Giant Stream - a stellar stream, with a well-defined 
eastern edge, which can be traced in projection to 
~ 70 kpc from the center of the galaxy towards the 
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Fig. 1. — A map of the surface density of blue RGB stars in a 
40n° area around M31 as mapped with the INT WFC. Overlaid are 
the locations of the HST/ACS point ings presented her e (red) and 
the deep HST/ACS halo pointing of IBrown et alJ 1^00^ ) (cyan). 



Radial velocity data have ruled out a simple connection 
to the dw ar f elli ptical companions M32 and NGC 205 
(p>ata et alH2004 (Fon t et al. 200l : 

NE Shelf - a large overdensity north-east of the galaxy 
centre which also exhibits a sharp outer boundary; 
NGC 205 Loop - a small (~ 15 kpc) loop of stars which 
appears, in projection, to emanate from NGC 205. 
While the geometry is very suggestive of a connection 
between the two, extant ste llar kinematical d ata still 
leave open other possibilities (|McConnachie et^a l. 2004); 
Gl Clump - a stellar overdensity which sits 30 kpc 
along the south-western major axis, in close proxim- 
ity to the anomolous globular cluster Gl. Although 
initially a connection between the two was sus- 
pected, radial velocity data have now ruled this out 
(iReitzel. Guhathakurta. RichI 120041: iFereruson et al.1 



Norther n Spur first hinted at by 

IWalterbos &: KennicuttI (|1988D . this is an extended 
plume of stars which sits just above the north-eastern 
major axis. 

Each field was observed for one orbit in the F606W 
filter (broad V) and two orbits in the F814W filter (I), ah 
within a single visit. Integer pixel dithers were executed 
between multiple sub-exposures to facilitate warm pixel 



and cosmic-ray rejection. Given the low stellar density 
in the outer halo, two separate pointings were obtained 
in each of the three outermost fields (Minor Axis, NE 
Shelf and Giant Stream) in order to ensure sufficient 
statistics on the upper RGB and horizontal branch. 

Images were first processed through the ACS pipeline 
and those in a g:iv en passband were combined using the 
Multidrizzle task (|Koekemoer, Fruchter, Hook, <fe Hac^ 
I2002D within PyRAF. We used the default mode of Mul- 
tidrizzle which calculates the offsets between dithered ex- 
posures using the information in the image header world 
coordinate system. The final drizzle was conducted using 
the LanczosS kernel with pixfrac and scale set to unity. 

Photometry was subsequently obtained usi ng th e 
IRAF implementation of DAOPHOT (Stetson 1987^. 
The stellar density in our fields is low enough to allow 
precise results from aperture photometry alone. PSF- 
fitting photometry was also carried out to compare with 
the aperture photometry and to measure how accurately 
individual sources could be modelled by a stellar PSF. 
For this purpose, a spatially non- varying PSF was built 
using ~ 150 bright stars in each of the nine ACS point- 
ings. Sources were retained in our final photometry list if 
their magnitude errors, sharpnesses and chi values all lay 
within 3-sigma of the average value at their magnitude. 
Aperture corrections were derived from several tens of 
bright stars in each combined frame and our photometry 
was placed on the STMAG system using the latest ze- 
ropoint vahies ([g irianni et al.Jl2005| h . Inspection of the 
iSchlegel. F inkbeiner. Davis' ('1998^ reddening map in- 
dicates values of E(B— V)= 0.06 — 0.08 (with uncertainty 
16%) in the far outskirts of M31; we adopt a constant 
value of E(B— V)= 0.07 for all fields and values of the 
selective extinct ion coefficient in the F606W and F814W 
passbands from Sirianni et m (j200 5). The final CMDs 
contain between 10,000 stars (the Minor Axis field) and 
90,000 stars (the Northern Spur field) and extend to 
four magnitudes below the horizontal branch (see Fig- 
ure 2). Artificial star tests indicate 75% completeness to 
F814Wo = 28.5 in fields of average stellar density. De- 
tails of our analysis will be presented in a forthcoming 
paper. 

3. RESULTS 

3.1. The CMD Morphologies 

Figure 2 shows Hess diagrams for the six halo 
fields targetted in our survey. In all cases, a broad 
RGB and a prominent red clump are evident, consis- 
tent with the fin dings of earlier HST studies of th e 
M31 halo (e.g. "Holl and. Fahlman. R,iche3 (|T99l : 
Bellazzini et al. (2003)). However, on closer inspec- 
tion, noticeable differences can be seen in the morphol- 
ogy of many of the CMDs. For example, the Giant 
Stream and NE Shelf exhibit rather prominent RGB 
bumps (F814Wo ^ 26.1). Note that the lower stel- 
lar density of the Minor Axis CMD makes it difficult 
to assess whether the RGB bump is similarily strong 
in this field. The RGB bump arises due to an evo- 
lutionary pause on the first ascent RGB when the H- 
burning shell passes through the discontinuity left by the 
deepest penetration of the convective envelope (Thoma^ 
119671: llbeni 119681) . It is a short-lived phase and, al- 
though commonly seen in single- age populations, has 
only been recently identified in composite populations 
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Fig. 2. — Hess diagrams of six halo fields in M31 shown with a square root stretch. Gl Clump, NGC205 Loop, Northern Spur, Giant 
Stream and NE Shelf are regions of stellar substructure within the M31 halo, while Minor Axis was chosen to represent clean halo. Clear 
differences are apparent between most CMDs, indicating genuine stellar population variations within the halo substructure. 



(e.g. iMonacQ. Ferraro. Bellazzini. Pancinol (|2QQ2f) ). 

The clarity of this feature in the shelf and stream CMDs 
suggests these regions contain a dominant stellar popula- 
tion with a small to moderate range of age and metallic- 
ity, perhaps the result of one or more bursts. Addition- 
ally, the location of the RGB bump at approximately 
0.4 magnitudes fainter than the red clump in these fields 
supports the inference from the RGB colour that this is a 
relatively metal-rich population (e.g. lAlves Saraiedinil 

The Gl Clump is distinct in showing clear evidence 
for a young population. It exhibits a well-populated 
upper main sequence (blue plume), with some stars as 
bright as F 814Wn ^ 24.5 . Comparison to the Padova 
isochrones (" Girardi 1 120041) indicates the brightest stars 
are consistent with a p opulation of age of ~ 2.5 x 10^ 
years (Faria et al. 2005). On the other side of the ma- 
jor axis, and at a slightly smaller radius, the Northern 
Spur shows no evidence for such a young component. 
Instead, this field shows a low- amplitude RGB bump 
at F814Wo ~ 26 and another bump at ~ 24.5. Given 
the moderate metallicity inferred from the RGB colour, 
the models of Alves & Saraiedini (1999) indicate that the 
brighter bump is more likely to be an asymptotic giant 
branch (AGB) bump as opposed to a secondary RGB 



bump (such features are also marginally present in the Gl 
Clump CMD). The NGC 205 Loop is similar to the Giant 
Stream, NE Shelf and Minor Axis fields in having an 
extended blue horizontal branch, but it differs in having 
neither an RGB /AGB bump nor a blue plume. Finally, 
all three CMDs in the top panel of Figure 2 show tan- 
talizing evidence of a bright main-sequence turn off at 
F814Wo ^ 27.5, roughly matched by a popula tion of 
- 2.5 X 10^ years with [Fe/H]= -0.7 (iGirardi II200I . 
As the Giant Stream and NE Shelf fields contain more 
stars than either the Gl Clump or NGC 205 Loop, the vis- 
ibility of this putative turnoff is not simply dependent on 
total stellar density. 

Although the CMD morphologies often differ, the 
mean colour of the (broad) RGB remains surprisingly 
constant. When measured in the region 27 < F814Wo < 
26, the mean RGB colour varies by less than 0.05 mag- 
nitudes from region to region. This is comparable to 
the difference expected in colour between [Fe/II]= —0.4 
and —0.7 stellar isochrones for a population of 14 Gyr 
(;Girardil[2004) . 



3.2. Using Stellar Populations to Trace the Giant 
Stream 
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The striking resemblance between the Giant Stream 
and NE Shelf CMDs motivates further exploration of the 
relationship between these features. Figure 3 shows the 
F814Wo luminosity function (LF) of these populations in 
the vicinity of the red clump and RGB bump. We fit the 
LF in the red clump region with the sum of a first order 
polynomial and a gaussian, selecting stars with the same 
range of de-reddened magnitude and colour in each case. 
We find the peak magnitude of the red clump to differ 
slightly between the fields, being F814Wo = 25.67 ±0.01 
for the stream and 25.53 ± 0.01 for the shelf. 

The magnitude and colour of the red clump are 
know n to vary as a f unction of age and metallicity 
(e.g. iQirardi <fe Salaris| (pOOl)) but the overah similar- 
ity of the CMD morphologies suggests that population- 
driven differences between the stream and shelf regions 
are small. This is further confirmed by the nearly iden- 
tical luminosity function slopes observed. In this case, 
a more likely cause of the offset in red clump magni- 
tude is that the regions are located at different line-of- 
sight distances. The red clump magnitude difference of 
A{F8UWo)sTREAM-SHELF = 0.14 ± 0.01 implies that 
the Giant Stream regi on lies beyond the N E Shelf by 
a factor of 1.07 ± 0.01. iMcConnachie et all ^003) have 
measured the line-of-sight distance to the stream as a 
function of position along it. We have interpolated their 
measurements to the location of our Giant Stream ACS 
field and find a line-of-sight distance of 830 ± 20 kpc. 
Combining this value with the differential distance mod- 
ulus inferred from the red clump magnitudes yields a 
distance to the shelf of 776 ± 20 kpc. 

4. DISCUSSION 

Our HST/ACS survey of M31 substructure has re- 
vealed significant variations in CMD morphology at large 
radius, indicating that the age and/or metallicity of stars 
in these parts is not constant. This is the first direct ev- 
idence for large-scale population inhomogeneities in the 
outskirts of M31. Prior HST studies, which have inferred 
a uniform population, have generally been based on ran- 
dom pointings within the inner h alo and hav e missed the 
main stellar overdensities (e.g. IBellazzini et al.l ((200 3)). 

It is tempting to interpret the observed population 
variations as evidence for multiple satellite accretions 
within the M31 halo, each having a distinct star forma- 
tion and chemical evolution history. One puzzle, how- 
ever, is the rather uniform mean RGB colour (though 
broad distribution) observed from field- to-field. If in- 
terpreted purely in terms of metallicity, the mean RGB 
colour suggests a variation of ^ 0.3 dex throughout the 
substructure. This could suggest that much of the halo 
debris comes from the accretion of a single system - 
one which has undergone little chemical evolution dur- 
ing several orbits around M31 but experienced at least 
one burst of star formation. On the other hand, age 
and metallicity may be varying in tandem throughout 
these regions, conspiring to produce a remarkably con- 
stant mean RGB colour (for examp le, as seen in the C a- 
rina dwarf spheroidal galaxy (e.g. iRizzi et 
We will evaluate these scenarios in more detail in subse- 
quent papers where we will model the CMDs in terms of 
star formation and chemical evolution histories. 

Previous studies of the M31 inner halo have demon- 
strated that ~25-50% of the population is old and 
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Fig. 3. — F814Wo LPs in the vicinity of the red clump and RGB 
bump for the Giant Stream (left) and NE Shelf (right). Overlaid 
(dashed line) are gaussian fits to the red clump. 



metal poor f e.g. IHolland. Fahlman. RicheJ (|1996D : 
Brown et al. (200^). The presence of a blue horizontal 
branch in four of the six CMDs in Figure 2 suggests this 
smooth halo population is also present at larger radii. 
The Northern Spur and Gl Clump CMDs, which probe 
the most extreme radii in our sample, do not, however, 
show this feature. Assuming our Minor Axis field is rep- 
resentative of the smooth halo component in M31, we can 
ask whether we should see this component in the outer- 
most fields studied. We find that, for both a spherical 
and flattened (b/a= 0.6) R^/^ distribution, the amount 
of smooth halo expected in the Northern Spur and Gl 
Clump CMDs is almost entirely negligible compared to 
the number of stars in the overdensity (4-10% of the de- 
tected stars). In other words, the absence of an obvious 
extended blue horizontal branch in these fields does not 
necessarily mean that the smooth halo component seen 
at smaller radii is not also present. 

The currently-favoured orbit for the stream progenitor 
is highly radial and passes very close to the g alaxy center 
before re-emerging nor th-east of the nucleus (|Ibata et al.l 
120041: iFont et al.l l2004l. Although radial velocities have 
not yet been published for stars in the NE Shelf re- 
gion, the strong similarity between the stream and shelf 
CMDs found here appears to s upport this predic tion. 
Furthermo r e, the calculations of llbata et al.l (|2004D and 
iFont et all (|200 4') indicate that the shelf region probed 
by our ACS pointing should lie many tens of kpc in front 
of our stream pointing. Our differential distance esti- 
mate based on the red clump magnitude is 50 ± 30 kpc 
a nd is thus in e ncouraging agreement. 

iBrown et al.l ( 2003) have recently obtained a very deep 
CMD of the M31 halo at a location 11 kpc along the 
southern minor axis (see Figure 1). The data reveal 
the surprising discovery of a significant (30% by mass) 
intermediate-age (6-8 Gyr) metal-rich population - a re- 
sult somewhat at odds with the traditional picture of 
stellar halos as old and metal-poor. The results pre- 
sented in this Letter suggest caution in interpreting this 
finding, as well as all others that are based on observa- 
tions of small single fields. The stellar populations in 
the outskirts of M31 exhibit large-scale variations and 
hence inferences about halo formation drawn from such 
pointings may not be representative of the entire halo. 
This is especially true in the case of Brown's field, which 
lies in close proximity to two regions which are strongly 
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contaminated by stream debris. 
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